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1. Introduction 

Matrix models have been extensively applied to various domains of physics and mathe- 
matics. Matrix integrals may primarily be viewed as tools for generating possibly decorated 
graphs of given topology (see also e.g. the reviews 0] and Q). The physical interpre- 
tation of these decorated graphs allows to view them as configurations of some statistical 
models with fluctuating "space" , namely defined on random lattices. This framework has 
provided a host of toy models for discretized quantum gravity. Indeed, physicists have 
engineered matrix integrals to suit their needs: the particular choice of integrand sets the 
peculiarities of each model, possibly including interactions between site variables living on 
the graphs. One then asks questions regarding thermodynamic and critical properties of 
large graphs, and much work has been done toward the identification of critical universality 
classes for the various models at hand. 

Graph combinatorics has been a subject per se ever since the ground-breaking work of 
Tutte, regarding the enumeration of various classes of planar graphs 0. The remarkable 
feature of most of these results is the existence of algebraic functional relations for the 
generating functions involved. When physicists started using matrix models to generate 
planar graphs, some of these results and relations were immediately identified through the 
algebraic character of the matrix model solutions. Somehow the simplicity of Tutte's results 
bore a "matrix model flavor". These methods always involved computing integrals over 
N x N Hermitian matrices in the large iV limit, and the main trick consisted in translating 
the set of integration rules into a diagrammatic technique, eventually leading to planar 
graphs. On a different front, some of these results have recently been revisited || in the 
light of an almost generic relation between this class of problems and that, much simpler, of 
enumeration of (possibly decorated) trees. By engineering suitable bijections, this allowed 
to reinterpret many of these functional relations in terms of trees, and to extend some of 



them so as to include extra colorings • It was first noticed by Morris |10[ that square 
complex matrix integrals could be used to generate so-called "chequered" surfaces, namely 
graphs with bicolorable faces, in such a way that no two faces sharing an edge be of the 
same color. More generally, we will show that the natural tools for enumerating face- 
colored graphs with prescribed numbers of faces of each color are integrals over rectangular 
complex matrices. 

Rectangular matrix models were first introduced in the context of nuclear physics 



by Bonk [ I I ] and first reemerged in an attempt to reach new critical models of two- 



dimensional quantum gravity p!]. More recently, rectangular matrix models were used 
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to calculate the partition function for effective Quantum Chromodynamics with arbitrary 



numbers of flavors and quark masses [13]. Rectangular matrix ensembles have also been 
studied extensively in the mathematics literature, where they bear the name of complex 
Wishart ensemble, or Laguerre ensemble ]TJ|] ||15|| , and recently found some combinatorial 



application to the study of increasing subsequences in random permutations [16 1 . The 



combinatorial application we have in mind in the present paper is of a very different 
nature. In Ref. ||12|| , the graphical rules for matrix integration were shown to generate 
polymer-like discrete surfaces. Here we will see that these rules naturally lead to discretized 
surfaces in the form of face-colored random graphs of fixed topology. Our aim is to apply 
known results from (rectangular) matrix integration to the actual computation of various 
counting functions involving face-colored graphs. 

The paper is organized as follows. In Sect. 2 we first recall the link between Hermi- 
tian matrix integrals and graph combinatorics, and extend it to the case of rectangular 
complex matrices in relation to colored-graph combinatorics. In Sect. 3, we recall a few 
facts on rectangular matrix integration, including the reduction to eigenvalues, various 
methods of solution and give a first combinatorial application of the Gaussian rectangular 
matrix integral to the enumeration of face-bicolored graphs of fixed genus with only one 
vertex. In Sect. 4, we present the computation of the general one-rectangular-matrix in- 
tegral. We first obtain the all-genus free energy using orthogonal polynomials, and then 
go on to the direct derivation of the planar limit, using a saddle-point technique. Sect. 5 
is devoted to the problem of enumeration of vertex-tricolored random triangulations via 
a rectangular multi-matrix integral. We show that the problem reduces to a single rect- 
angular matrix integral with a logarithmic potential, and recover the result first found 
in [17] using a two-Hermitian matrix model and then rederived in a purely combinato- 



rial manner in 0. Section 6 gathers a few concluding remarks and two possible further 
applications of rectangular matrix integrals to statistical models on random tessellations, 
namely Interaction-Round-a-Face (IRF) models and hard-particle models. 



2. Matrix Model Combinatorics 

As a preliminary example!, let us consider the N x N one Hermitian matrix integral 

JdMe- NT ^-^>-^ , 
2? = i (2 1) 

JdMe- NT ^ 

2 We refer the reader to the lecture notes || for a more pedagogical and thorough introduction 
to the subject of matrix model combinatorics. 
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Fig. 1: Pictorial representation of Gaussian Hermitian matrix integration. 
The matrix elements are represented by oriented double half-edges (a), 
whose lines carry the matrix indices i,j respectively. A typical vertex repre- 
senting Tr(M 4 ) is represented in (b): matrix indices are conserved along the 
oriented lines. The result of the Gaussian integration is to connect all the 
half-edges into pairs forming edges (c) along the lines of which the matrix 
indices are conserved. 




Fig. 2: A typical planar graph obtained from the one Hermitian matrix 
model integral. The pairs of lines bordering the edges are oriented and carry 
matrix indices, which must be summed over from 1 to N, resulting in a weight 
N per face of the graph. 

where the measure of integration is the standard flat Haar measure dM = Tii<j dMijdMij Y\ i dM„ 
for which the denominator of ( |2.1|) plays the role of normalization factor. This ratio of 
integrals may be computed by formally expanding the g^-terms of the exponential, and 
by explicitly performing the term-by-term Gaussian integration. The Gaussian integrals 
are alternatively calculated by the following graphical procedure. To each matrix element 
Mij one associates an oriented double half-edge as shown in Fig. 1 (a), where each line 
carries an index of the matrix entry. In this representation, a factor of the form Tr(M fc ) is 
pictured as a /c-valent vertex with k outgoing half-double edges as displayed in Fig. 1 (b) 
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for k = 4; these vertices carry the weight Ng^. By virtue of Wick's theorem, the result of 
the Gaussian integration is obtained by summing over all possible ways of connecting the 
half-edges of the integrand into pairs so as to form a closed graph; the corresponding prop- 
agator (gluing factor) (MijMki) = Sudjkt/N is displayed in Fig. 1 (c), and is accompanied 
with a weight inverse of the prefactor of the Gaussian term in (|2.1|) , namely t/N. Note 
that on the resulting graph (see Fig. 2 for an illustration) the oriented lines forming the 
edges now carry an index i — 1, 2, N of the matrix, which must be summed over in the 
end, resulting in a weight iV per loop of the graph, i.e. equivalently a weight N per face 
of the graph. The computation of (|2.1|) therefore reduces to drawing all possible fatgraphs 
with oriented double-lined edges like that of Fig. 2, and counting them with a weight Ngk 
per /c-valent vertex, k = 1,2,..., as well as t/N per edge and N per face, together with 
an overall inverse symmetry factor coming from the various combinatorial numbers when 
expanding the exponential and resumming contributions of equivalent graphs. Moreover, 
taking the logarithm of Z restricts this sum to only connected graphs, receiving therefore 
an overall contribution N 2 ~ 2h , where h is the genus of the graph. Hence expanding the 
free energy Log Z at large N yields a formal power series generating arbitrary graphs with 
fixed topology and valences. The graph of Fig. 2 clearly has h = and contributes to the 
planar free energy. 

The reduction of the integral (|2.1] ) to one over eigenvalues, together with standard 
orthogonal polynomial techniques allow for a complete although tedious calculation for all 
N. For combinatorial purposes however it is interesting to concentrate on planar (genus 
zero) graphs, selected by taking the limit N — > oo. In this limit, the integral ( |2.1| ) can be 
computed at leading order by use of a saddle-point approximation, and allows in particular 
to recover the generating function for rooted Eulerian maps (graphs of even valences, 
namely #2i-i = and g2i ^ in (|2.1|), and with a marked edge), first obtained by Tutte 
|J. The computation of the planar limit of ( |2.1|) in the case (fe-i = for all i is detailed in 
appendix A for completeness. Note that a complete combinatorial picture involving trees 



was recently found in [18] for the general case of arbitrary g^s, including odd indices. 

In this note we will concentrate on the combinatorial applications of rectangular ma- 
trix integrals. For the purpose of illustration, we now consider an integral over complex 
rectangular matrices A of size N\ x N2 of the form 



(p© 



A 



A — ^r® i v 

®<S) 

(a) (b) (c) 

Fig. 3: Pictorial representation for Gaussian rectangular matrix integra- 
tion. The matrix elements Aij,A\- are represented (a) as double half-edges 
with oriented index lines marked 1,2 for the corresponding range of indices 
Ni, N2. An extra overall orientation distinguishes between A and A* . In (b) 
we have represented a typical vertex for Tr(AA^) 2 : the oriented lines carry 
the corresponding conserved indices. The result of the Gaussian integration 
is to connect the half-edges A to the At into pairs through edges of the form 
(c); the indices are conserved along the oriented lines. 

where the measure of integration is flat over the complex elements: dA = Yli j dAijdAij. 
Note that the integrand is unitary invariant like in ( |2.1| ): this is what makes the integral 
easy to calculate, by changing variables to the eigenvalues of A At. We compute Z again by 
expanding the g^-terms in the exponential, and by performing their term-by-term Gaussian 
integration. The latter is again obtained through a similar graphical procedure to that of 
the Hermitian matrix case. We have represented in Fig. 3 respectively (a) the double 
half-edge for a matrix element A^-, carrying the indices i = 1,2, Ni and j = 1, 2, N2, 
where the distinction between the two ranges of indices can be made by coloring the two 
corresponding oriented lines say with colors 1 and 2; (b) the 2/c-valent vertex corresponding 
to a term Tr((AAt) fc ) for k = 2 and (c) the propagator (A^A^) = 5u5jkt/N. Using again 
Wick's theorem, each term in the expansion is integrated by connecting all half-edges into 
pairs so as to form a closed graph like that of Fig. 4. Note that the form of the propagator 
allows only to connect lines with indices running over the same set. Summing over those 
results finally in a weight iVj per corresponding loop of color i. Equivalently, this amounts 
to coloring the faces of the graph with colors 1 and 2 (white and grey in the graph of Fig. 
4), and to weight them with N±, N2 respectively. Writing Ni = Nxi for convenience, and 
expanding Log Z, we get the generating function for arbitrary face-bicolored graphs with 
a weight t per edge, X{ per face of color % = 1,2, gk per 2/c-valent vertex, k = 1, 2, an 
overall factor of JSf 2 ~ 2h ^ fi the genus of the graph and finally the usual inverse symmetry 
factor. 
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Fig. 4: A typical planar graph obtained from rectangular matrix integration. 
The edges of the graph are made of two oriented lines, each carrying a matrix 
index running respectively from 1 to Ni and N2. The faces of the graph are 
colored in white (resp. grey) according to the index range N\ (resp. N2) of 
the oriented line bordering it. With the above rules, face colors must alternate 
around each vertex. Summing over running indices results in a weight Ni per 
white face and N2 per grey one, namely a weight iVf iVj for the present graph. 

It is now clear how rectangular matrices will help solve colored-graph combinatorial 
problems. But before we go into this, we need a bit more experience with rectangular 
matrix integration, the subject of next section. 



3. Rectangular Matrix Models 



In this section we recall a few known facts \T2] on complex rectangular matrix inte- 
gration, and apply them to the warm-up exercise of counting face-bicolored graphs with 
one vertex. 

3.1. Reduction to eigenvalues 

As opposed the case of Hermitian (hence square) matrices, where the first step of 
integration consists in reducing the integral over the matrix to one over its eigenvalues, 
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this is no longer possible with a rectangular matrix A as such. However both AA^ and 
A^A are Hermitian with non- negative eigenvalues, and one may still reduce the rectangular 
matrix integrals at hand to integrals over these non- negative eigenvalues. 

More precisely, let us consider the general unitary invariant rectangular one-matrix 
integral 



Z N {V-N U N 2 ) 



dAe -NTrV(AA*) 



(3.1) 



where the integral extends over the complex rectangular matrices A of size N\ x N2, and 
we use the flat measure dA = Yli=i Y[j=i^i,j^i,j- We have picked the normalization 
of the measure arbitrarily, as we will only be interested in computing ratios of objects of 
the form ( |3.1|) corresponding to different Vs. Let us assume iVj < AT 2 . Diagonalizing the 
Ni x Ni Hermitian matrix AA^ leads to AA^ = PA 2 P^, where A = diag(Ai, AjvJ, a 
diagonal matrix with non-negative entries, and where P G U(Ni) is unitary. Repeating 
this for A' A leads to another diagonal matrix containing the A 2 's and completed by N2—N1 
zero eigenvalues. Upon rearranging the eigenvalues by permutation, we may write A* A = 
Q(A') 2 Q^ where A' = diag(Ai, A^, 0, 0) and Q G U(N2). Upon picking positive 
square roots of the eigenvalues and modifying P and Q accordingly, we finally arrive at 
the "rectangular diagonalization" formula 



A = PXQ\ X 





V 





Ao 













A AT, 







0/ 



(3.2) 



where X has size N\ x N 2 . It is then tempting to change the variables of integration in 
( pTT| ) to A — > (P,Q,A), and noting that the integrand only depends on A, to integrate 
out the angular variables P, Q. This can indeed be done but we must be careful with the 
change of variables itself, which is not bijective as such. Indeed, if $ denotes an arbitrary 
diagonal matrix of size Ni x N± with diagonal entries of the form e Wi , and R an arbitrary 
unitary matrix of size (A^ — N%) x (A^ — N±), we see that the transformations 



P -> P$, 



$ 
R 



(3.3) 



leave A unchanged. Therefore the right locus for the angular integration is (P, Q) G 
(U(N X ) x U(N 2 )/U(N 2 - Nx)) /U(l) Nl . The Jacobian of the transformation is then easily 
computed (see the appendix of Ref. |I(J for the case Ai = N 2 , and Ref. for the 
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straightforward generalization to N± > N 2 ) and after integrating out the angular variables 
we get up to some proportionality constant (the volume of the angular variables locus) 

Z N {V;N U N 2 ) = f J] jf d\ t \^e- NV ^A A(A) 2 (3.4) 

where A (A) = rii<i<,/<jVi (A? ~~ * s the Vandermonde determinant of A. The two 
notable differences with the Hermitian case are that (i) the Aj's are all non-negative and 
(ii) the potential receives an extra logarithmic piece from the Jacobian, namely W(X) = 
V(A) — N2 ~^ Nl Log A, and this is the only place where N 2 appears, while Ni is the number of 
integration variables A^. It is worth noting that the form of the potential W obtained here 
is similar to that of the so-called "generalized Penner model" studied in Ref. |TI| (with 
a parameter t = — (N 2 — Nx)/N) where however it is used to define (square) Hermitian 
matrix integrals. 

3.2. Methods for solving 

The next step uses the well known property of the Vandermonde determinant that for 
any momc polynomials p m (A) = A m +0(A m " 1 ) we have A (A) = det(A^ 1 ) = det(p 3 -_i(Ai)). 
We may now conveniently pick the p m 's to be orthogonal with respect to to the one- 
dimensional measure d\w(X) = d\\ N2 ~ Nl e~ NV ( x \ namely 

n OO 

{p n ,Pm)= / d\\ N *- N ^e- NV ^p n (\)p m (\) = h n 5 n , m (3.5) 
Jo 

for some normalization factors h n , and finally expanding the determinants one gets 

jVi-l 

Z N {V]N 1 ,N 2 ) = N 1 l J] hj (3.6) 

3=0 

The partition function (V; iVj , N%) is then explicitly computed by deriving recursion 
relations between the h^s from their definition ( |3.5| ). In the following sections, we will 
solve some explicit models for combinatorial purposes. 

In this paper we will be focusing on combinatorial applications of the partition func- 



tion (|3.1| ) for various choices of the potential V, by interpreting them as suggested above 
as generating functions for colored graphs. In particular, we will be interested in the limit 
of planar (genus zero) graphs, obtained in the large limit as explained above. Intro- 
ducing Xi = Ni/N, for i = 1,2, the planar limit is obtained by sending N —>■ 00, while 
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xi,X2,V remain fixed, and we will assume that x\ < x 2 . By a slight abuse of notation, 
we denote by Zn(V;xi,x 2 ) = Zn(V; Ni, N 2 ). We whish to compute the planar free en- 
ergy for the correctly normalized model, namely with partition function ( |2.2| ) expressed as 
Zn(V;xi,x 2 )/Zn(Vo;xi,x 2 ), where Vq(x) = x/t is the pure Gaussian potential. Defining 
the following planar free energies: 

/ = f(V;x 1 ,x 2 ) = lim — ^LogZ N (V;x 1 ,x 2 ) ^ 
N^oo iv z (3.7) 

fo = f(V ;x ll x 2 ) 

the correct generating function for bi-colored planar graphs is simply / — Jo- 
in the large iV limit, Zn(V;xi,x 2 ) may also be computed by a saddle-point method. 
Indeed it reads Zn(V;xi,x 2 ) = J d\i...d\N 1 e~ N s '( a i.--,^jv 1 ) j w here the action functional 
S is defined as 

1 Nl f 1 1 

S(X U X Nl ) = - J2\ v ( x >) ~ («2 - ^i)Log Xi I - —2 L °SI^ " ( 3 - 8 ) 

i=l ^ ' l<i^j<iVi 

The large N behavior is then dominated by the solution to the saddle-point equation 
dS/dXi = 0, namely 



Let us as usual introduce the resolvent 



1 Nl I 

1=1 



expressed in terms of the solution to ( |3.9|) . When expanded at large z, this function 
generates the quantities 

u ~ 1 (Tr((AA^)) - 1 ! dAT?{{AA^r)e- NV ^ 

u Ntm = -{Tr((AA) ))-- J dAe -NV(AAt) t 3 - 11 ) 

through un(z) = J2 m >o u N,m z ~ m i with leading term un,o = Ni/N = X\. Equivalently 
we may define the density of eigenvalues Pn(X) = -h Yl!i=i ^(^ ~Xi), a distribution related 
to the resolvent through 

Pn{X) = ^-(u N (X + iO) - u N (X - i0)) (3.12) 
9 



Multiplying both sides of ( |3.9| ) with l/(z — Xi) and summing over i, we arrive at the 
following quadratic differential equation for the resolvent: 

1 cIujn(z) , 2 



N dz 



+ u N (z) 2 - W'(z)u N (z) + P N {z) = (3.13) 



where W(z) = V (z) — (a; 2 — £i)Log z is the total potential including the extra logarithmic 
term from the Jacobian, and 



Vl W'{z)-W'{X i ) 



z — Xi 

( 3 - 14 ) 

= — Un{0) + z-X t 

When iV — > 00, the derivative term in ( |3.13|) vanishes and we end up with a quadratic 
equation for the resolvent uj(z) = liniTv^oo ujn(z), easily solved as 

2 

where the ± sign is selected by demanding that co(z) ~ x±/z at large z (a property clearly 
inherited from the definition (|3.10| )). To further fix co(z) completely, we will make the 
standard "one-cut" hypothesis, that this complex function only has one cut in the complex 
plane, translating into the fact that the limiting density of eigenvalues p(X) = limpN(X) 
has a compact support made of a single real interval. 

To finally recover the planar free energy (|3.71 ), let us note that / = — limjv^oo S(Xi, X^), 



evaluated on the solution of the saddle-point equation (|3.9|) , we then easily compute 

df dS v-k OS dXi 1 , 1 ^1 , 

t-r = -t-z- - > T^t—^- = -lim— > Xi = — 3.16 

dt dt dXi dt t N ^ t y ! 

i=l 1 i=l 

where we denote by iv m = limjv-^oo UN >rn , uJN,m as in ( ft.ll| ). So finally the corresponding 
generating function for bi-colored graphs with a marked edge reads 

E(t) = t f t ~ t ^ = \^ V ) - "iWO) (3-17) 

where the mention of V and Vq simply refer to the potential one must use to compute the 
corresponding coefficient of z~ 2 in the resolvent's expansion at large z. 

In the following sections, we will complete the calculations outlined here in the case 
of some specific potentials V. Let us for starters deal with the Gaussian case V — Vq. 
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3.3. A first combinatorial application: bicolored trees and the Gaussian model 

The simplest rectangular matrix model is the Gaussian one, say with potential 

Vq(AA^) = AA^/t (3.18) 



The corresponding orthogonal polynomials p n (satisfying (|3.5| ) with V = Vq) are expressed 
in terms of Laguerre polynomials and lead to a trivial result for Zn(Vq; x±, x%) oc 
is however instructive to derive first the large N (planar) limit of the model, to see in par- 
ticular how the Wigner's semi-circle limiting distribution of eigenvalues for the Hermitian 
case is modified in the complex rectangular case. Moreover, the corresponding resolvent is 
the generating function of face-bicolored planar graphs with one vertex, in bijection with 
bicolored trees. 

The quadratic equation ( |3.13| ) for the resolvent reads in the large N limit 



u 2 - (- - + P(z) = (3.19) 

Moreover we have 

, Nt 1 

P N (z) = X -^ 1 i- T V i- = -^—^ w ^(0) - P(z) = -^-J^ U {0) (3.20) 

TV An 
1=1 



so that finally 




X2 — X\ Ifl X2 — Xi\ 2 X2 — X\ 



, + 4^-^(0) (3.21) 

z V V z / z 

The negative branch of the square root has been selected to ensure that u(z) ~ x\/z for 
large z, which further fixes uj(0) = —l/(t(x2 — x±)) and finally 

u(z) = — (z- t(x 2 - xt) - ^z 2 - 2t(xi + x 2 )z + t 2 (x 2 - xi) 2 ) (3.22) 

The corresponding limiting density of eigenvalues reads p(A) = (uj(\+i0) —uj(\—i0)) / (2in) , 
so that 

p(A) = -^—^2t{ Xl + x 2 )\ - A 2 - t 2 {x 2 - xi) 2 (3.23) 
ZttiA 

This distribution has a compact support A G [(y/tx2 — \Ztx~i) 2 , (y / tx2 + y/txi) 2 ], away from 
the origin as X2 > X\. The corresponding planar free energy fo obeys eqn. ( |3.16| ), with 
the result 

= ^ = x lX2 (3.24) 
dt t 
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Note that in the limit x^ — > x\ — > 1, upon performing the change of variables A = y 2 the 
measure p(X)d\ reduces to the Wigner-semicircular distribution, with support [— 2\/i, 2y/t], 
The large N resolvent u(z) = J2k>o UkZ ~ k ~ 1 ma y be interpreted as the generating 
function for some bicolored graphs as follows. Indeed, in the planar limit, the coefficients 
iVk = lim7v^ 00 (Tr(M fc ))/A^ are obtained by drawing all planar, face-bicolored graphs, with 
one 2/c-valent vertex, and with a weight X{ per face of color i, and t per edge. Expanding 



u>(z) we recover the result of [|17 



K3 C™ O^^ 1 (3 - 25) 

m=0 \ / \ / 

where these graphs were identified as bicolored arch configurations, equivalent by duality 
to bicolored rooted trees. 

3.4- All-genus Gaussian model: face-bicolored graphs with one vertex 
We now wish to compute the all genus expansion 

k - — {Tr{AA ) ) = — NTrAAt/t 1^ N ° k > h ^ 3 " 26 ) 

J h>0 

where 9k,h is the number of genus h face-bicolored graphs with one 2/c-valent vertex, and 
with a weight t per edge and X{ per face of color i = 1, 2. In the previous section, we have 
computed the large iV limit u>k = 9^,0 ( |3.25|) . To go beyond this requires to go deeper into 



the orthogonal polynomial solution of the model. The orthogonal polynomials p n for the 
Gaussian case are nothing but a rescaled version of the generalized Laguerre polynomials 
L"(x), orthonormal with respect to the measure dxx a e~ x on IR_|_. More precisely, we have 
the generating function ||20|| : 



G( ,, ,) 3 my ,» (3.27) 

n=0 ^ ' 

with p n as in fl3.5|) , with V = Vq: the orthogonality property is readily obtained by 
computing 



poo ~n„..m / M\ n ^ m 

/ dxx N ^e- N -/ t G(x,z)G(x,w)= V ^ (p n ,p m ) 

J ^-^ n\m\ \ t J 

n,m>0 ^_ 28 ^ 

t/N 



n,m>0 
N 2 -N 1 +l 

1 — 2iy 
12 



The fact that the right hand side of ( p. 28 ) is a function of zw implies the orthogonality, 
and this immediately yields the normalization factors h n upon expanding the right hand 
side of (|3.28| ) with the result 

i , x N 2 -N 1 +2n+l 

hn=ljj) nl^-TVi+n)! (3.29) 

which finally yields 

Z n (V ; Xl , x 2 )= Nl l J] hn=U) m\ n n\{N{x 2 - x x ) + n)\ (3.30) 

n=0 ^ 7 n=0 

where we recover the expected t NlN2 behavior. 

To compute 9k, we use the following generating function 

H(u) = £ Vyp.0 h = i(Tr( e^V*)) ( 3 .31) 

fc>0 

Reducing to eigenvalues, it is a simple exercise to show that for any function /, the average 
(/) = (Si=a fi^i)) 1S expressed in terms of the orthogonal polynomials as 

(/) = e' < 3 - 32 > 

Applying this to /(x) = e Nux ^, we are left with the computation of the modified norms 
h n = (e Nux / t p ni Pn), while h n = (p n ,Pn) are given by (|3.29| ). Using again the generating 
function G(x,z) ( |3.27| ), we write 

[dxx N ^e N ^/ t G(x,z)G(x,w) = T Z ^(e N ^p n ,p m ) 
I n\m\ \ t I 

J n,m>0 7 



(N 2 -N 1 )\ 



(3.33) 

t/N xv, v,,. 



1 — zw — w(l + z)(l + 
Picking the diagonal terms in the z, w expansion, we finally get a closed formula 

V 7 m>0 r=0 V 7 V y 

which together with ( |3.29| ) and ( ^.32| ) allows to compute the generating function H(u) of 
( CT) as: 

n=0 0<r<m<r+n \ / \ / \ / 
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and finally picking the coefficient of u k , we get 



9 k 



k\t k 



tE e 



n=l r=max(0,k+n — Ni) 



JVi - 71 

k — r 



N 2 — n + r 
r 



(3.36) 



A first remark is in order. If we suppress the coloring, by taking Ni = N 2 = N, we get 
a different result from that of Ref. [^TJ obtained in the context of modular group com- 
binatorics, enumerating the all-genus graphs with one 2/c-valent vertex: indeed, arbitrary 
graphs with one 2/c-valent vertex are not necessarily bicolorable as soon as their genus 
h> 1. 




Fig. 5: The only genus one face-bicolored graph with one hexavalent vertex 
has one face of each color, and corresponds to the term X1X2/N 2 in #3. 



The first few numbers 9b read 



01 = xix 2 

6> 2 = XiX 2 (xi + x 2 ) 

9 3 = xix 2 {x\ + 3xix 2 + x 2 + "779 ) 

5 

4 = x\x 2 {x\ + x 2 )(xf + 5xix 2 + x\ + — j) 

5 8 

05 = x±x 2 (x 1 + 10xfx 2 + 20xfx 2 + 10xix 2 + x\ + -77^(3x1 + 8X1X2 + 3xo) + — -j) 

N 2 iv 4 

(3.37) 

As an illustration, we have represented the first genus one contribution appearing in 03 at 
order 1/N 2 in Fig. 5. 
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4. The one-matrix model: vertex-bicolored graphs 

Let us now consider the case of a general potential 



M 



V(AA*) = V^(AA t ) fc (4.1) 

k=i 

The combinatorial interpretation of Zn(V; xi, x 2 )/Zn(Vq\ x±, x 2 ) is clear: it generates ar- 
bitrary face-bicolored graphs with weights gk per 2/c-valent vertex, k = 1,2, ...,M, t per 
edge and Xi per face of color i = 1,2. Taking a finite upper bound M is just an artifact to 
make calculations easier, but our result will extend to any formal power series involving 
all g^, k > 1. The section is organized as follows. We first solve the problem by use of 
orthogonal polynomial techniques. The result is a set of recursion relations determining 
the free energy of the model, that generates connected face-bicolored graphs of arbitrary 
genus. In the subsequent section, we concentrate on the planar (large N) limit and derive 
closed expressions for the corresponding generating functions. 

4-1. Orthogonal polynomial solution 

We wish to compute the all-genus free energy 

F N (y;x 1 ,x 2 ) = -r^Log 



N 2 \Z N (V ;x 1 ,x 2 ) , 
= J2N- 2h f h (V;x 1 ,x 2 ) 

h>0 

where fh is the generating function for connected face-bicolored graphs of genus h with a 
weight t per edge, gk per 2/c-valent vertex and Xi per face of color i = 1,2. The relation 
( [OP allows to compute completely in terms of the normalizations h n (V),h n (Vo) of 



the corresponding orthogonal polynomials (3J3). To determine relations between the h's, 
we compute 

ot at t Z /. o\ 

d dp n N ^ ' 

Qj_(Pn,Pn-l) = (-Qj-,Pn-l) + ^ [xp n ,Pn-l) = 

The operator of multiplication by x acting on the p n 's is self-adjoint and therefore reads 

XPn(x) = p n +l(x) + S n p n (x) + r n p n _i(x), T n = (4.4) 
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Moreover, if we write p n (x) = x n — A n x n_1 + 0(x n ~ 2 ) and substitute it into ([4.4]), we 
easily identify X n = Y^j=o s j- Finally, by using the orthogonality of the p's, we see that 
(^§f-,p n ) = and that (^r,p n _i) = —\ n h n -\. This allows to rewrite ( fOj ) as 



^-Logh n = s n 

N 2 dt X n-2^ N 2 dt S l- r " 
3=0 



(4.5) 



Introducing the "partial" free energy F n = (1/N 2 ) J2o<j< n -i Log/ij, we may write r n = 
(t 2 d t ) 2 F n /N 2 which finally leads to a closed equation for F n , by noting that F n+ i +F n _i — 
2F n = (Logr n )/iV 2 : 

Together with the initial conditions F$ = and i 7 ! = (l/iV 2 )Log ho, this determines _F n 
completely as a power series of t, and allows to compute the desired all- genus free energy 
F N (V;xi,x 2 ) = (F n (V) - F n (y ))\ n=Nl , by using eqn. ( ^6|) in both cases of V and Vq. 
More precisely, let us introduce the subtracted partial free energy f n = F n (V) — F n (Vo). 
Using the last equation of (|4.5|) for both V and Vq, we may write 



r n (V) (t 2 d t ) 2 F n (V) 1 , (t 2 d t ) 2 f n 



Pn = 7TTT = x9 ^ , T , x = 1 + 



r n (V ) (t 2 d t ) 2 F n (V ) (t 2 d t ) 2 F n (V ) (4.7) 

= e JV 2 (/„ + i+/„_!-2/„) 

We now use eqn. ( |3.29|) to explicitly compute the partial Gaussian free energy F n (Vo) = 
(l/N^^JhoghM) as 

n-l 

F n (V ) = (1/N 2 ) Log {(t/N) N ^+ 2l+1 t\(N 2 -N 1 + (4.8) 

i=0 



and finally the denominator in (4~7) reads 



= (4.9) 
which leads finally to the following recursion relation for f n : 

f u/ •>/ 1 T , N 2 d t t 2 d t f n \ 
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to be supplemented by the following initial conditions 



/o = 



1 / [™ dxx N *- N *e- Nx /*e Nl Ei 9iS £\ (4.H) 
h = —Log ' I 



N 2 M f™dxx N »- N ie- 



■Nx/t 



Note that the details of the potential, in particular the dependence on the g^s is entirely 
contained in the initial condition Q4.11| ), which is supposed to be formally expanded as a 
power series of the g^s. 



To solve (|4.10|) ( [OTP , let us proceed order by order in t. Indeed writing 



1 oo 



N 2 
k=i 



then eqn. ( |4. 1U| ) implies the following recursion relation, upon expanding the logarithm: 



f n+ i, k + U-i, k - 2f n , k = J2 J X n)p E UWi + Wn,* (4-13) 

p>l F V 1 1 ' fe i >l,i = l,2,...,p i=l 



to be fed with the initial conditions (|4.11|) that /o,fc = and that is identified with 



the coefficient of t k in the expansion of N 2 fi of ( f4.11| ). Moreover, it is easy to prove by 
induction that all f n ^ are polynomials of n of degree k + 1, and always have a factor 
n(N 2 - iVi +n). 

17 



Explicitly, for the first few values of k, we find 
N 2 f n ,! = n(N 2 -Ni + n) 9l 

N 2 f n , 2 = \n{N 2 -N l+ n) (g 2 + N * - + % 2 ) 
N 2 f n , 3 = \n{N 2 -N x +n) (gf + ^(N 2 - N t + 2n) 9l g 2 

+ ^ ((N 2 - iVi) 2 + 5n(iV 2 - N x ) + 5n 2 + l)g 3 ) 
N 2 f nA = \n{N 2 -N x +n) (g{ + ^(N 2 - N ± + 2n)g 2 g 2 

+ ^ ((N 2 - iVO 2 + 5n(iV 2 - N 1 ) + 5n 2 + l)g m 

+ ^ ((2(iV 2 - N x ) 2 + 9n(N 2 - N x ) + 9n 2 + l)g 2 

+ j^(N 2 -N l + 2n)((N 2 - N x ) 2 + 7n(N 2 - N,) + 7n 2 + 5)g 4 ) (4.14) 

N 2 f n , 5 = ^n(N 2 -N x + n)(gl + ^(N 2 - N, + 2n)gfg 2 
10 

+ ^ ((N 2 - iVO 2 + 5n(iV 2 - N x ) + hn 2 + l)g 2 g 3 

+ ^(N 2 -N 1 + 2n)((N 2 - iVi) 2 + 7n(N 2 - Ni) + 7n 2 + 5) 9l g 4 

+ ^ ((2(iV 2 - iVi) 2 + 9n(iV 2 - iVi) + 9n 2 + l) 9l g 2 

+ ^3 (N 2 - iVi + 2n)((N 2 - N x ) 2 + 6n(N 2 - + 6n 2 + S)g 2 g 3 

+ ^4 ((N 2 - N x ) 4 + Un(N 2 - N x f + (56n 2 + 15)(iV 2 - N x ) 2 
+ Un(6n 2 + 5)(N 2 - iVi) + 2(21n 4 + 35n 2 + 4))g 5 ) 





Fig. 6: Two genus one graphs corresponding to terms of order 1/iV 2 in -Fjv,4, 
with respective weights g\x\x 2 (2 four-valent vertices, one face of each color) 
and g±x\x\ (one eight-valent vertex, one white and two grey faces). 
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To get the final result for the subtracted free energy Fn(V; xi, x%) = J2k>i ^N,kt h °f 
I, we simply note that Fjv,fc = fn,k\n=Nn hence we finally get 

Fn,i = xix 2 gi 

F N ,2 = ^xix 2 (gf + (xi + x 2 )g 2 ) 

Fn,3 = 7^xix 2 (gf + 3(xi + x 2 )gig 2 + {x\ + 3xix 2 + x\ + -^Osfc) 

F N ,4 = ~XiX2 {gf + 6(xi + x 2 )g\g 2 + A(x\ + 3xix 2 + x\ + ^-) gi g 3 

+ (2x1 + 5xix 2 + 2x1 + tttt)^! + (xi + x 2 )(x\ + 5xix 2 + x\ + -77, W 

F N ,5 = ^x x x 2 (gl + 10(xi + x 2 )g\g 2 + 10(x 2 + 3xix 2 + x\ + j^)glgz 

+ 5(xi + x 2 )(x\ + 5xix 2 + xl + j-p)gig4 + 5(2x1 + 5x x x 2 + 2x\ + —)g ig * 

3 

+ 5(xi + x 2 )(x\ + Axix 2 + x\ + —pr)g 2 g3 + (x\ + 10xfx 2 

JM z 

+ 2§x\x\ + lOxixl + x| + ^(3x 2 + 8xix 2 + 3x1) + J^)9s) 

(4.15) 

The quantities Fjy,k are the generating functions for arbitrary face-bicolored graphs with 
k edges, and with a weight N~ 2h , h the genus, a weight Xi per face of color i, and a weight 
gi per 2i-valent vertex, and also an inverse symmetry factor. We have displayed in Fig. 
6 for illustration some of the genus one graphs contributing to i*j\r,4- When N — > oo, we 
may identify easily several coefficients in i*jv,fe- In particular, the coefficient of g\ simply 
counts graphs made of one loop with k marked points, i.e. with a symmetry factor 1/k, 
and therefore two faces, one of each color. Moreover, the coefficient of gk in Fn^ (all-genus 
graphs with one 2/c-valent vertex) is easily identified with 9k of ( |3.37|) . 

In the following section, we use an alternative method to directly get the planar limit 
of the free energy in a more compact form. 

4-2. Planar limit 

Repeating the computation of Sect. 3. 2, we finally arrive at the planar limit of the 
resolvent in the form 

. zV'(z) - x 2 + Xi + ^(zV'(z) -x 2 + Xi) 2 - 4z 2 P(z) 
u{z) = — (4.16) 

where zV'(z) — x 2 + xi is a polynomial of degree M in z, and zP(z) a polynomial of degree 
M — 1. Applying the "one-cut" hypothesis, let us assume that the polynomial under the 
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square root has M — 1 double zeros, and two single ones, say b > a > 0, so that the density 
of eigenvalues has the compact support [a, b]. This amounts to writing 



iV'(z) -x 2 + X! + Q(z)y/(z -a)(z- b) 
2z 



(4.17) 



where Q is a polynomial of degree M — 1. Imposing that u(z) ~ x\j z at large z, we may 
solve for Q: 

zV'(z) — X2 + Xi 



Q(z) 



(4.18) 



y/(z- a)(z - b) 

where the bracket must be expanded at large z, and the subscript + simply means that only 
the polynomial part of this expansion must be retained. Note that, with these definitions, 
the term — x 2 + x\ in (ft-lEp may be dropped as it only participates to terms of order 1/z 
and higher. By analogy with the Gaussian result ( |3.22j ), let us introduce the following 
parametrization for a and b 



a=(v^2- y/U^) : 



b = (Vu~ 2 + Vu~ 1 y 



(4.19) 



so that a + b = 2(Ui + U2) and ab = (U2 — Ui) 2 . Introducing the polynomial coefficients 
il>k(Ui, U 2 ) through the following large z expansion 



1 



m±M + Ml! ,Vo z 



m=0 ^ ' 



(4.20) 



and substituting V'(z) = 1/t — X1i<A;<m 9k z i n to ( [4.18| ), we get 



Q(z) 



Viz) 



1 a-\-b I ab 

z z 2 _i _|_ 

M-l M-l-m 



(4.21) 



-J+J2 Zm 9m+k+l^k{U U U 2 ) 



m=0 k=0 

To further fix a and 6, we must moreover write that oj{z) ~ X\jz at large z. Introducing 



(4.22) 



the large z series R(z) = —V'(z)/ J 1 — + ||, we may rewrite 

u(z) = (xi - x 2 + (Q(z) - i2(z))V(*-a)(z-&))/(2z) 
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so that the above condition reads x\ — x 2 — R(z)\-i = 2x\ (the subscript simply means 
that we must pick the coefficient of the corresponding power of z in the Laurent series for 
R), and eventually 

1 M 



(4.23) 



fe=i 



We must finally ensure that the resolvent uj is well-defined at z = 0, which amounts to 
imposing that the numerator of u(z) in ( |4.17| ) vanishes! at z = 0, i.e. that —VabQ(0) = 
X2 — X\, namely, using \fab = U 2 — U%: 



M-l 



X 2 



- Xl = (JJ 2 - U x ) 9k+iMUi,U 2 )) 



(4.24) 



k=0 



The two equations (|4. 23|) ( ^j2J ) determine V\ and U 2 (i.e. a and b) completely upon 
requiring that U\ , U 2 be formally expanded in powers of t around the origin, starting with 
Ui = txi + 0(t 2 ), i = 1, 2. These two equations may be recast into 

U x M 



k=l 
M 



x 2 

<Pk(Ui, U 2 ) 

where the (fk's are generated by 

1 1 + 



(4.25) 



fc-i 

E 

m=0 



k=l 



k\ fk- 1 



m 



m 



U1-U2 



2 I 1 _ 2 u 1 ±u z+ (Ui-u 2 r 



Efk_ 
yk 

k>0 



(4.26) 



Let us now express the free energy of the model as in ( |3.16|) t-4? = We must use 
again the series R(z) to express the 1/z 2 coefficient of u(z) in ( |4.22| ) as u\ = (—R(z)\- 2 + 
^R{z)\-x)/2. Together with (fyj|), this finally yields 

M 



dt 



t x t 



k=l 



m=0 V 1 / V / 



(4.27) 



3 We must pay attention here to the fact that ( 4.17 ) defines u correctly as a power series of z 
at infinity. Going to the limit z — > implies an overall change of sign of the square root term. 
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where the fa's are generated by 



2U X U 2 




U1+U2 



U1+U2 , (Ui 



-c/ 2 ) s 



E <Pk 



(4.28) 



k>\ 



We finally arrive at the complete generating function for rooted (edge-marked) bicol- 
ored graphs with a weight t per edge, Xi per vertex of color i = 1, 2 and g k per 2/c-valent 
vertex, k 



1,2,.., M: 



E(t) = ^-t 



4f d/ Z7iZ7 2 



A/ 



(4.29) 



where U\,Ui are the unique solutions to ( [4.25| ), such that U,, 
becomes even simpler if we compute the derivative 



k=i 

tXi + 0{f). The result 



d 



(tE(t)) 



U X U 2 



xix 2 



(4.30) 



dt"~ v "" t 2 

by using ( [125| ) to eliminate both dlli/dt and dUildt. 

A few remarks are in order. Eqns. ( |4.25| ) were first obtained in § in a purely 
combinatorial way, by establishing bijections between face-bicolored planar graphs and 
suitably colored trees, generated by the functions U\,U<i (denoted tP, tQ there, while xi, £2 
are denoted u, v). Although neither expression ( |4.29j ) nor ( |4.3(J| ) appear there, we checked 
their compatibility with the equation dE/dx2 = {U\/t) — X\ obtained there (actually 
expressed in terms of a function B = E + X1X2). Eqn. (|4.30|) can be rephrased as a 
bijection between suitably marked face-bicolored graphs and pairs of such trees. 

As a non-trivial check of our calculations, let us suppress the colors by taking x\ = 
x 2 = 1. We then have U x = U 2 = U, and i^ k {U, U) = C^ 1 )^ from ([0§). This reduces 
eqn. ( fl.25|) to eqn. (A. 13) of appendix A upon redefining M — > J and gi — > g2i- Noting 
that (f) k (U, U) = (j^)^" 1 we also find that (ggg ) reduces to eqn. (A.14) of appendix A. 
Finally, eqn. ( |4.3Q|) clearly reduces to eqn. (A. 15) of appendix A. 

As a final example, the generating function for rooted bicolored quadrangulations, 
namely with gj~ = 8^2, reads 



E(t) 



U1U2 
t 2 



(1 - 2t(Ux + U 2 )) - x x x % 



1 



xi = C^i (--(^i + 2t^)) 
X2 = U 2 {j-(U 2 + 2U 1 )) 



(4.31) 
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a,: 
a;: 



-0 

-6) 



4- 




(D© 



(a) 



(b) 



(c) 



Fig. 7: The pictorial representation for the triple Gaussian integration over 
Ai,A2, As with respective sizes N% x N2, N2 x iV~3, iV~3 x N±. Each matrix ele- 
ment is represented as a double half-edge (a), each oriented line carrying one 
index (we have indicated the color % = 1, 2, 3 of each index, corresponding to 
its range Ni, N2, N 3 respectively. Each half-edge carries an additional overall 
orientation, to distinguish between Ai and A\. We have represented the two 

vertices Th:{AiA2A 3 ) and Tr^gA^!) in (b). The result of the integration 
is to connect the half-edges into pairs by forming edges (c) along the lines of 
which indices (and therefore index colors) are conserved. 

5. More matrices: vertex-tricolored triangulations 

5.1. Rectangular multi-matrix model 



In ||17|| , the generating function for vertex-tricolored random triangulations with a 
weight t per edge of color (12), and Xi per vertex of color i = 1, 2, 3 was obtained by use of 
a two-Hermitian matrix model. Here we present an alternative solution, using the following 
integral over three rectangular matrices A\, A2, As respectively with sizes Ni x N2, N2 x N 3 , 

iv 3 x Nr. 



Z N (V;N l ,N 2 ,Ns) = J dA^dAse'^^^^^ 

V(Ax,A 2 , As) = -(AiA\ + A 2 A\ + A 3 A\) - A 1 A 2 A 3 - A\A\A\ 



The free energy of the model reads 

1 / Zm(V- iVi A^o Ni) 
F N (V; N lt N 2 , Ns) - F N (V ; N lt N 2 , N 3 ) = -^Log N ) / -± 2 - 3 

N 2 \Z N (y ;Ni,N2,N 3 ) 



(5.1) 



(5.2) 



where V (A 1: A 2 , A 3 ) = ^{AiA\+ A2A\ + AsA\) is the Gaussian part of the potential. The 
expansion of this free energy now involves summing over all connected trivalent fat graphs, 
with double-edges carrying indices ranging up to Ni, i = 1,2,3 according to the corre- 
sponding matrix indices. Indeed, expanding the non-Gaussian part of the exponential, 
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we must again integrate term-by-term. This is done diagrammatically using the elements 
of Fig. 7: (a) the matrix elements (A^jk are represented as double half-edges as be- 
fore, but now with three possible line colors 1,2,3 corresponding to indices running up 
to Ni, N 2 , N3 respectively; (b) the terms Tr(AiA 2 A 3 ) and Tr^gAjjjAj) correspond to two 
trivalent vertices around which colors cyclically alternate respectively in clockwise and 
counterclockwise order; (c) finally the half-edges around the vertices must be paired into 
edges through the propagators ((A m )ij(A m y kl ) = SuSj^r/N, in such a way that line colors 
are conserved. 

We therefore must draw all connected face-colored trivalent fatgraphs using colors 
1, 2, 3, such that the three faces adjacent to each vertex have distinct colors. These graphs 
are nothing but the duals of vertex-tricolored triangulations. Further setting X{ = Ni/N, 
i = 1, 2, 3, we see that these graphs receive a weight x\ per face of color % = 1, 2, 3, r per 
edge, and an overall jV 2-2 ' 1 , h the genus of the graph. The precise contact with the results 
of |17| is made through setting t = r 3 , as there are three times as many edges as those 
colored (12) in the triangulations. 

The integral (|5.1|) may be simplified by explicitly integrating out say A 3 . Indeed, 
forming the product 

-A 3 Al - A 3 A X A 2 - A\a[a\ = -(A 3 - tA\a[)(A\ - tA x A 2 ) - tA\a\A x A 2 (5.3) 
r ' r 

and performing the Gaussian integration over A' 3 = A3 — rA^Aj, we are left with a two- 
matrix integral with potential 

V(A 1: A 2 ) = -{AtAl + A 2 A\) - tA\A x A 2 A\ (5.4) 
r 

Now the integration over A 2 can also be explicitly performed, as a particular case of the 
general formula 

JdAe- T ^ AAH ^ _ 1 

JdAe- T ^ AAt ) ~ det(M) Af 3 (5 - 5) 
where A is rectangular complex of size ./V2 x ./V3 and M non-singular complex of size N 2 x N 2 . 
We are finally left with a rectangular one- matrix integral over A = rA\ of size N\ x N 2 , 



and we may simply substitute in the logarithm of (5^) the ratio of partition functions 
with Zn(V)/Zn(Vq) where 



Z N (V) = J dAe- NV ^ 



A /it 

V(AA^) = —— + x 3 Log(J - AA^) (5.6) 
Vo(AA^) = *f 
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where we have set t = r 3 . So we finally see from simple direct matrix integrations that the 
counting of vertex-tri colored triangulations reduces to that of bicolored Eulerian graphs, 
with a weight £3 per vertex of any (even) degree, namely taking infinitely many g^s in 
(fPD, and taking them all equal to x$. This relation was actually implicitly used in fl? 



where a two-Hermitian matrix model actually counted the relevant Eulerian graphs. A 
more general relation was found in [0 connecting the so-called /c-constellations to k + 1- 
valent Eulerian graphs, and the present situation corresponds to k = 3. 

5.2. Orthogonal polynomial solution 

At this point we just have to use the results of Sect.4.1, by simply setting g^ = X3 
for all k > 1 in the various equations i.e. with V(x) = x/t + X3Log(l — x). Some of them 
however must be handled with care. Indeed, the main recursion relation ( |4.10|) remains 
the same, as it basically is independent of the particular form of the potential V, but the 
initial condition ( f4.11| ) is changed. It now reads 

/o = 

1 fj^dxx^-^il-xr^e-^/^ 
fl ~ iV^ L ° g ^ dxx^-Nie-N*/* ) (5.7) 

/ OO / ~\T , 7 1 \ /AT" ~\T 1 7 \ / , \ 

^ Log E 



N 3 + k — 1\ (N 2 - N i + k \ k] (_L 



iV 2 \^ V k J V k J \N 

\k=0 v 7 v 7 v 

where the term (1 — x)~ Na is only meant as a formal power series of x. Again, the final 
result for the free energy is obtained by taking f n \ n= N 1 - 

It is quite interesting to compare this result for the all-genus free energy of the vertex- 



tricolored triangulations to that of [pT7[| , reading with the present notations 



b n+1 + n _i - 2</> n = ^Log (l + —(td t ) 2 c 



N 2 ' 




■n 



(5.J 



1 t ^fN 3 + k-l\fN 2 + k-l\,.ft 



and the final result for the free energy is obtained by taking 4> n \ n =N 1 as before. That ( |4.10| ) 
and (|5.7| ) lead to the same result as (|5.8|) , namely that 4 > n 1 = /jVd is a highly non-trivial 
fact. 
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5.3. Planar limit 

The planar limit of the free energy is easily obtained from the results of Sect. 4. 2, again 
by first having infinitely many g^s in the potential (|4.1|) , and by finally taking = x 3 for 
all k > 1. Using the generating functions (|4.26| ) (|4.28|) of the cp^s and </>fc's at z = 1, and 
substituting the result into the extended versions of Ql29| ) and ( J4.25| ), we get 

^ = ^ (j - f ( = - 1) ) (5-9) 







-u 2 






- ?7i - J7 2 


?- 

u 2 


AU X U 2 




j/i - ^ 
i+c/2 - 


Ui 






U x - U 2 ) 


2 - 


AU X U 2 



We then use the following change of variables to get rid of the square roots: setting 
Ux = Vi(l - V 2 ) and U 2 = V 2 (l - V x ) indeed transforms (1 - U x - U 2 f - AU X U 2 into a 
square (1 — V\ — V 2 ) 2 . This results in 

E(t) = ^((1 - Vi)(l - V 2 ) - 1 _y'_y 2 ) ~ X lX2 

tx 1 = V 1 (l-V 2 - p ) (5.10) 
1 — V\ — v 2 

Introducing finally a third variable V3 = £3:3/(1 — V\ — V 2 ), we finally recover the result of 



E(t) = ^^(l-V 1 -V 2 -V 3 



tx 2 = V 2 (l-V 1 -V 3 ) 
tx 3 = V 3 (l-V 1 -V 2 ) 

where Vi = txi + (9(£ 2 ) are to be expanded in power series of £ and substituted into the 
expression for E. This result was given a purely combinatorial proof in Ref. ||, by suitably 
mapping marked tricolored triangulations to colored trees generated by the V^'s. 



6. Rectangular multi-matrix models for new random lattice statistical models 

In this paper, we have shown how complex rectangular one-matrix models could be 
used to generate and count various types of face-colored graphs. The methods of solution 
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presented are adapted from the Hermitian matrix case, and lead in particular to closed 
formulas for the planar generating functions, expressed in terms of tree-like generating 
functions, thus providing an alternative proof of some recent combinatorial results based 
on some bijections between rooted colored planar graphs and some suitably colored trees. 

The standard methods of Hermitian matrix models can also be adapted to the rectan- 
gular multi-matrix case, by use of the Harish-Chandra-Itzykson-Zuber formula |[22|1 , namely 
that 



r , ts det \e aih A^. . ... 

/ dUe^ aUbU ^ oc L A/ J?f:^ N (6.1) 

Ju(N) A(a A(6) 



'U(N) 

where the integral is performed over unitary N x N matrices and a, b are some diagonal 
matrices of same size. An interesting generalization of this formula was also derived in 
Refs. P3] pi p5[ for rectangular matrices and reads: 



I dill dVe-^-'^c **V»r*Wl l <i J <s, (6 . 2) 



where A and B are rectangular matrices of size N\ x N2 (with N\ < N2), a and b are the 
positive square roots of the diagonalized versions of AA^ and B^B respectively, and I n are 
the modified Bessel functions. Rectangular multi-matrix integrals with chain-like interac- 
tions are expected to generate more involved face-colored graphs, and can be reduced to 
eigenvalue integrals by use of either ( |6.1| ) or ( |6.2| ) depending on the form of the interaction 
terms (AA^BB^ or AB* + BA*). 

We give below two series of examples of new rectangular multi-matrix models, whether 
solvable or not, describing random lattice statistical models. 

6.1. IRF Models on random tessellations 

Hermitian matrix models are well adapted for generating configurations of vertex 
models, as the terms in the potential are directly related to vertex configurations. Here 
we concentrate on so-called Interaction-Round-a-Face (IRF) statistical models, originally 
defined on fixed lattices by assigning height variables to the vertices of the lattice (say a 
map i — ► h(i) from vertices to a target space T), and attaching to each configuration of 
these heights a Boltzmann weight, itself a product of local Boltzmann weights depending 
only on the heights around each elementary face of the lattice. A tremendous amount 
of work has been done for solving fixed lattice IRF models, in relation with integrability. 
When defined on a random lattice, these require that the faces of the random tessellation 
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replacing the lattice keep their original form, so that the local Boltzmann weights remain 
the same. A nice global picture for a particular class of random lattice IRF models (the 
so-called Restricted Solid On Solid models based on the A,D,E Dynkin diagrams) was 
obtained in PBfl , using first the correspondence between these models and loop gases in the 



planar case, and then using a method of surgery to generate higher genus contributions. 






8, 



8, 



g, 



Fig. 8: The four vertices of the matrix model describing the general two-state 
IRF model on a random triangulation. These vertices are dual to the faces of 
the triangulation. According to the nature of the matrix indices running over 
the lines (the small circled numbers 1, 2 stand for the ranges j = 1, 2, N± 
or N2 respectively), the faces (i.e. the vertices of the dual triangulation) 
may be colored accordingly with 1,2 (the large circled numbers 1,2). The 
corresponding Boltzmann weights #1,2,3,4 are indicated. 



Let us show how rectangular matrix models naturally generate general IRF models. 
The trick is to use a dual picture of IRF models, in which the height variable is a face 
variable (say a color) defined on the irregular faces of a random lattice with regular vertices. 
The particular choice of potential fixes the corresponding vertex Boltzmann weights, which 
dually translate into the local face boltzmann weights of the IRF model. For simplicity, 
we start with IRF models on random triangulations. We use a rectangular matrix model 
to generate the dual (trivalent) graphs, for which the original height variable h(i) simply 
becomes the color of the face z, itself indicated by the range of matrix indices running 
along the loop bordering the face j = 1, 2, iVj. For instance, in the case of a target with 
two elements, say T = {1,2}, we must assign colors 1,2 to the faces of random trivalent 
graphs, with specific Boltzmann weights #1,2,3,4 depending on the three colors around each 
vertex, as shown in Fig. 8. This is readily realized by a triple integral over two Hermitian 
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and one rectangular matrices A of size Ni x N\, B of size N 2 x N 2 and C of size N\ x A^: 

_ J dAdBdCe~ NTr ( v ( A ' B,c ^ 
Z ~ J dAdBe~ NTl ( v °( A > B > c » 

V (A,B,C) = ^- + ^ + C& (M 
F (A, B, C) = V (A, B,C)- 9l ^--g 2 ^-- g 3 AC& - g 4 BC^C 

Note that the potentials V and Vq only make sense here within a trace, as A and S have 



different sizes. As before, the logarithm of Z in (|6.3|) generates the configurations of the 
most general 2-state IRF model on triangulations of arbitrary genus h, with a weight 
A 7 " 2-2 , provided we have set A^ = ccjiV. Particular cases of this model include: (i) the 
triangular Ising model, with N\ = e 11 N, N 2 = e~ H N, gi = g 2 = ge 3K / 2 , gs = g$ = 
ge~ K l 2 1 where H is the magnetic field, K = J/(ksT) the normalized spin coupling and 
g the cosmological constant and (ii) Baxters's triangular three-spin model with an 
additional magnetic field H, with Ni = e H N, N 2 = e~ H N, g\ = g 4 = ge K , g 2 = g% = 
ge~ K . As the matrix C occurs only in quadratic terms in the potential V of (|6 . 3|) , it 
may be integrated explicitly Indeed, writing this quadratic form as Ylijki ^ij^iki^jk^U ~ 
g^SjkAu — g^BjkSu), the integration over C allows to rewrite the model as a two-Hermitian 
matrix model, with matrices A and B of different sizes N\ and A^, and with potentials 

A 2 R 2 
V (A,B)= T + T 

1 A 3 B 3 

V(A, B) = V (A, B) + —Log (I Nl ® I N2 - g 3 A ® I Na - g 4 I Nl ® B) - 9l — - g 2 — 

(6.4) 

where the subscript of the identity refers to the size of the matrix, and again the potentials 
only make sense within a trace as A and B have different sizes. This extends easily to 
the case of /c-state IRF models. Note that in the particular case of A, D, E-type targets, 
where states of adjacent vertices of the triangulation may either be of the same color 
or of neighboring ones on a vertex-colored target graph of A, D, E-type, the rectangular 
matrices can all be integrated out, leaving us with a slight generalization of Ref. [27] in 



which matrices may have different sizes, which has the effect of introducing a "magnetic 



field" coupled to colors. Let us also mention the square version of (|6.3| ) for completeness. 
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It generates the most general two-state IRF model on random quadrangulations: 

_ f dAdBdCe- NT *( v ^ B > c » 

Z ~ J dAdBe- NT *( v °( A > B > c » 

A 2 B 2 
V (A,B,C) = — + — +Ctf 

V(A, B, C) = V (A, B,C)- 9l ^-g 2 ^- g 3 A 2 C& 



(6.5) 



2 



- 94 B 2 tfC - g 5 ACB& - ge^p- 

Note that if g$ 7^ 0, the rectangular matrix C can no longer be integrated out as in the 
previous case. 

It is now clear that increasing the number of matrices will allow for exploring larger 
target spaces T. One other case of interest is the straightforward generalization of the 
vertex-tricolored triangulations of Sect. 5 into random /c-gonal tessellations with vertices 
colored 1, 2, k cyclically around each face, clockwise or counterclockwise in alternance, 
weighted by £Cj per vertex of color % and t per edge. The planar free energy for this 
problem was derived recently in and is also related to work on constellations J7|. The 
k colors require the introduction of k rectangular matrices A\, A2, A^ of respective 
sizes Ni x N2, N2 x ./V3, Nk x N\, while the coloring rule is encoded in two (clockwise 
and counterclockwise) vertices A\Ai---Ak and its adjoint. The partition function for this 
problem reads 



Z 



fnLi dA * e ~ NTi(Vo{Ai '-' Ak)) 

A t A\ (6.6) 



V = V - g(AiA 2 ...A k + 44_ 1 ...A t i ; 

Taking again Ni = XiN yields the desired all-genus free energy F = LogZ. The particu- 
larly simple form of the planar free energy found in suggests that this model should be 
solvable. 

6.2. Hard objects on random tessellations 

Another class of interesting models is that of mutually excluding (so-called hard) 
particles. On a fixed lattice, vertices may be empty or occupied by a particle, with the 
constraint that no two particles may occupy adjacent vertices, and one moreover attaches 
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an activity z per particle. The case of the triangular lattice was solved by Baxter [BS 
under the name of hard hexagons (the area excluded by a particle on the triangular lattice 
is the hexagon formed by its six adjacent vertices). Neither the square nor the hexagonal 
lattice cases were solved exactly yet although much is known on them [[29] []30 



in par- 



ticular thanks to the powerful corner transfer matrix method developed by Baxter pi 



As far as critical properties are concerned, all these models are expected to undergo a 
crystallization transition at some positive value z = z c of the activity, between a fluid low- 
density phase and a crystalline high-density phase, in which particles tend to maximally 
occupy a sublattice of the lattice. Remarkably, this clearly distinguishes the triangular 
lattice case from the two others, as there are three candidate sublattices for the crystal 
as opposed to only two in the square and hexagonal cases. As a consequence, the hard 
hexagon model lies at the crystallization transition in the universality class of the critical 
three-state Potts model (unitary CFT with central charge c = 4/5), while in the two other 
cases, the expected universality class is that of the critical Ising model (unitary CFT with 
central charge c = 1/2). 
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Fig. 9: The three vertices used via matrix model to generate the dual graphs 
to the configurations of the hard particle model on random quadrangulations. 
We have indicated the nature of the matrix index lines by a small circle (white 
for indices ranging from 1 to N, grey for indices ranging from 1 to P = zN). 
We have also represented in thin lines the squares dual to the vertices, and 
the corresponding status of their (empty or occupied) vertices by large (white 
or grey) circles. The Boltzmann weights <7i 2,3 f° r the three situations are also 
indicated. 



Recently, it was shown in [|32[ that the expected critical properties for square and 
hexagonal lattices are correctly reproduced on random vertex-bicolorable lattices, which 
still allow for the existence of two crystalline states. These involved however lattices with 
faces of arbitrary valences. Rectangular matrices can be used to generate analogous models, 
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Fig. 10: A typical configuration of the hard-particle model on the faces of 
a random tetravalent planar graph, as obtained from some particular mixed 
square and rectangular matrix integral. Occupied faces are represented in 
grey, while empty ones are white. These correspond respectively to indices 
ranging from 1 to P and N on the oriented line bordering them. When 
summed over, these result in a weight N per empty face and P per occupied 
one. By duality, this graph is nothing but a random quadrangulation with 
hard-particles at some of its vertices. 

but with regular faces, the randomness being concentrated into the arbitrary valences 
of vertices. For the case of hard squares, we must generate random quadrangulations 
with empty and occupied vertices according to the exclusion rule. Note that the vertices 
are naturally bicolored, and that we therefore expect to reproduce the correct critical 
behavior of hard squares, as we will always have two natural maximally occupied crystalline 
states. This requires simply two matrices, one square Hermitian A of size N x N and one 
rectangular B of size N x P, whose integral will generate the duals of the decorated 
quadrangulations. Indeed, using the three vertices displayed in Fig. 9, we may attach 
to the faces of the tetravalent graph they form once connected two colors say white and 
grey (empty or occupied) according to whether the running matrix index around the face 
ranges from 1 to N or P, as illustrated in Fig. 10. The allowed vertices of Fig. 9 carry 
the information on the exclusion rule. The partition function for the hard particle model 
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on random quadr angulations then reads: 

_ f dAdBe- NTl W A > B » 
Z ~ j dAdBe- NT < v ^ A - B )) 



V Q (A 1 B) = ^-+BB^ ( 6 - 7 ) 
V(A, B) = V (A, B)- 91 —- g 2 A 2 BB^ - 9s [ -^- 



Indeed, upon setting P = Nz, the logarithm of Z of (|6.7|) generates the all-genus connected 
quadrangulations with hard particles on their vertices, and a weight _/V 2_2/l , h the genus, a 
weight z per particle and extra weights g\ per empty square, g 2 per singly occupied square 
and gs per doubly occupied square (which can be interpreted as a diagonal interaction), 
and finally the usual inverse symmetry factor. The integral (|6.7|) may be transformed into 
one over two Hermitian matrices as follows. We first form the square 

-g 2 A 2 BB^ - g 3 { -^- = -%(BB* + —A 2 ) 2 + ^-A± (6.8) 
2 2 # 3 2g 3 

which we represent by the use of a iV x N Hermitian matrix M as 

e _^(-^( BB t + || A 2 )2) = JdMe 

JdMe~ NT ^ 

Integrating over B is readily done using ( |5.5|) , and now yields an additional term zLog(l — 
yTg^M) in the potential. We therefore end up with an integral over the two Hermitian 
matrices A, M of size N x N, with potentials 

2 Jl,f2 



A 2 M 



2 2 



2X A 4 (6-10) 



V = V - -^=MA 2 - ( 91 -2^)4- + ^Log(l - v^M) 

x/93 V 03/4 



This is a particular generalization of the 0(n = 1) model coupled to discrete 2D quantum 
gravity. When gig 3 = 2g 2 , it indeed reduces to a particular 0(n = 1) model, with a 
logarithmic potential. We expect its (multi-) critical point to be that of the critical Ising 
model coupled to 2D quantum gravity. 

We may now also write a candidate model for the random surface version of hard 
hexagons, by considering the hard particle gas on a vertex-tricolored triangulation. Indeed, 



we expect just like in [32] that preserving the tricolorability of the triangulation will allow 
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Fig. 11: The allowed vertices for the random lattice version of the hard 
hexagon model, via rectangular matrix integral. The small white (resp. grey) 
circled numbers i = 1,2,3 refer to matrix indices ranging from 1 to Ni (resp. 
Pi = zNi), while the large white (resp. grey) ones refer to the color of the 
empty (resp. occupied) corresponding face or vertex of the dual triangula- 
tion. We have indicated the nine different rectangular matrices needed for 
the model. 

for the existence of three competing crystalline states, eventually leading to the critical 
transition of the three-state Potts model on a random surface. The idea is simply to 
attach a different color to occupied colored faces, say by replacing the range of matrix 
indices j = 1,2,..., Ni (for an empty face of color i) by j = 1,2,..., Pi = zNi (for an 
occupied face of color i). The exclusion rule is then enforced through the allowed vertex 
configurations of Fig. 11. The model reads as follows: 



Z 



I nti dA l dB l dC l e- N ^ v ^> B M» 



V = Y,A i A\ + B i B i i +C i C, 



(6.11) 



V = V - g(A 1 A 2 A 3 + B 1 A 2 C 3 + C 1 B 2 A 3 + A t C 2 B 3 + h.c.) 
where h.c. denotes the hermitian conjugate of the sum of four terms (the total of eight 
terms corresponds to the vertices of Fig. 11), and where the nine rectangular matrices 
have respective sizes 

A x : Ni x N 2 A 2 : N 2 x N 3 A 3 : N 3 x N\ 

B 1 : zN x x N 2 B 2 : zN 2 x N 3 B 3 : zN 3 x JVi (6.12) 
d : iVi x zN 2 C 2 : N 2 x zN 3 C 3 : N 3 x zN x 
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The logarithm of Z (|6.11|) generates the configurations of the hard-particle model on 
vertex-tricolored random triangulations of genus h, with a weight J\f 2 ~ 2h provided we take 
Ni = XiN, i = 1, 2, 3, and also with a weight z per particle and x^ per face of color i. The 
partition function Z of (|6.11|) may be further simplified by explicitly integrating out (i) 
the Q, i = 1,2, 3, which results in new terms -g 2 {BiA 2 At,B\ + B 2 A 3 A\b\ + B^AXbI) 
for the potential and finally (ii) the Bi, % = 1, 2, 3, by use of the formula (|5.5| ). The final 
result is an integral over the three rectangular matrices A^. 



j dA l dA 2 dA^e- NT ^ vi ~ A ^ A ^ 
~ J dA 1 dA 2 dA 3 e- NTr ( v °(AuA2,A3)) 

V = A X A\ + A 2 A\ + A 3 A\ 

V = V - g(A 1 A 2 A 3 + 444) + ^iLog(l - g 2 A 2 A\) 
+ 2x 2 Log(l - g 2 A 3 Al) + zx 3 Log(l - g 2 A 1 A\) 



(6.13) 



This is easily generalized to the hard-particle model on the vertices of a cyclically 
/c-colored random /c-gonal tessellation (with alternating clockwise and counterclockwise 
order), with the additional constraint that a given /c-gon have its vertices occupied by a 
total of at most one particle. Such a model naturally gives rise to k distinct crystalline 
states. The corresponding rectangular multi-matrix model is reduced in an analogous way 
to an integral over k matrices Ai, A 2 ,... A^ with sizes N\ x N 2 , N 2 x N 3 , ... Nk x N±, and 
with potentials 

k 

i=l 

k s k \ 

V = V Q - g(A l A 2 ...A k + 4^!-i— 4i) + Yl ^ Lo §( 1 ~9 2 ( H Aj) x h.c. ) 

3jii— l,i— 2 

(6.14) 

where we have extended the definition of indices in such a way that = k and — 1 = k — 1 
(namely when % = 1 the terms A k -i and A k must be omitted in the product, while when 
i = 2 A\ and A k must be omitted). In the case k = 4 of hard particles on cyclically vertex- 
4-colored quadrangulations with at most one particle per square, we expect the model 
to reach a crystallization transition in the universality class of the random lattice critical 
4-state Potts model (CFT with central charge c = 1). Note that this point is not reachable 
within the model ( |6.7|) of hard-particles on plain quadrangulations with g 3 = (i.e. with 
at most one particle per square), for which no natural maximally occupied crystalline state 
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will survive. More precisely, taking g% = in ( |6.7|) allows to integrate over the rectangular 
matrix B, leaving us with a one-Hermitian matrix model with potentials 

V (A) = ^, V(A) = ^- + ^Log(l - g 2 A 2 ) - 9l ^- (6.15) 

and it is a simple exercise to check that the only multicritical points of this model all 
correspond to negative critical values of z and therefore lie in the universality class of the 
random lattice Lee- Yang edge singularity (non-unitary CFT with c = —22/5), meaning 
that the (unitary) crystallization transition has been wiped out by the sum over random 
quadr angulations. 
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Appendix A. Planar Eulerian Graphs from the one-Hermitian Matrix Model 

We wish to compute the large iV limit of the function Z (|2.1|), in the even potential 



case, namely with V(x) = X/i=ifl'2i%^- The cutoff degree 2 J is purely artificial, but 
dealing with polynomials greatly simplifies the calculations. As we are only interested in 
formal power series, we will simply take J —>■ oo in the final result. To compute Z of 
( |2.1|) , we first reduce the matrix integral to an eigenvalue integral by changing variables 
M — > (U,m), where M = UmU' , where m =diag(mi, m_/v) and U G U(N)/U(1) N 
are some angular variables. The Jacobian of the transformation is well-known to be the 
Vandermonde determinant squared A(m) 2 , with A(m) = Yli<j( m j ~ m i)- One is therefore 
left with the computation of 

. Z(V) 



^ Z(V ) 
Z(V) = J dme- N2s ^ 



1 N I 
= —J2 V ( m i) ~ J^2 Yl LogK-m,, ^ v 

i=l l<i^j<N 



V (x) 



x 
2t 
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We now evaluate separately the numerator and denominator in the large N limit, where 
the corresponding integrals are dominated by the saddle-point of the action functional S, 
namely the m^'s such that dS/drrii = 0, hence 

V'(m i ) = ^- y y (A.2) 

v lJ N ^ m i -m j y ! 

This leads to a quadratic equation for the planar resolvent 

1 - 1 

u(z) = lim — Y (A.3) 

i=i 

evaluated on the solution m to ( |A.2| ), and related to the limiting density of eigenvalues 

N 

;Y — -x /V 



1 - 

P(z) = (A.4) 



i=i 



through p(z) = (u)(z — iO) — uj(z + i0))/(2i7r). The solution takes the form 

u(z) = \{V\z) + ^V'(z) 2 -AP(z)) (A.5) 

for say g2i > 0, where P(z) = limjv^oo jt ^2 i= i(V' '(z) — V (mi)) / \z — mi) is a polymomial 
of degree M — 2. lo(z) is further fixed by its large z behavior u(z) ~ 1/z and by also 
demanding that it have only a one-cut square root singularity in the complex plane. This 
amounts to saying that V — 4P has M — 2 double zeros and two single ones say a, 6, so 



that finally uj(z) = (V'(z) + Q(z)^(z — a)(z — 6))/2, Q a polynomial of degree M — 2. 
Moreover, the potential V(x) being an even function of x, it is clear that the density of 
eigenvalues p is also an even function, and therefore b = —a. The (even) polynomial Q is 
fixed by requiring that the large z expansion of oj(z) have no positive powers of z. This 
gives 



Q{z) 



V'(z) 



(A.6) 

+ 



\J z 2 — a 2 

where we must expand the bracket at large z, and the subscript + indicates that one only 
retains the polynomial part of the expansion. Finally, the value of a is determined as 
follows. Writing R{z) = —V'(z)/\/z 2 — a 2 as a series expansion at z = oo, we also get 



that u)(z) = — \ \R{z)\ _ \/ z 2 — a 2 , where the subscript — means that we only retain the 
negative powers of z in the expansion. As V is odd, the constant term z° vanishes and we 
must simply express that the residue (1/z term) of u(z) is equal to 1, giving the relation 
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R-2 = —2, where the subscript indicates that we only retain the corresponding power of 
z term in the expansion of R. Explicitly writing 



J z 2 - a 2 ^ z 2k + x z 2k + x 

v * " fc>0 fc>0 v 7 



and further defining Ik = for k < 0, we get 



1 T 



fc=i-j i=i 



which, applied to the relation i?_ 2 = —2, finally leads to the equation determining a: 



The single zero a of the eigenvalue density is the unique formal series solution to (|A.9| ) 
that starts as a = 2\/t(l + 0(t, {<?2i})) for small positive t and g's. 

Finally the planar free energy reads f(V) = lim7v^oo(l/-^ r2 )LogZ'(y) = — S{m, V) 
taken at the saddle-point solution. We also have 



df dS drrii dS , 1 v-^ <^2 

t-*- = -y- — t— L -t—= lim — = — A.10) 

dt ^ dun dt dt n^oo N ^ 2t 2t K ' 
i=i i=i 

where we have expanded the resolvent as u>(z) = 1/z + z 2 + c^/z 3 + 0(l/z 4 ) at large 
z. We may therefore identify tdf /dt with 

t |=I(-2 fl _ 4 -2^) (A.11) 

For the Gaussian potential Vo, the equation ( |A.9|) reduces to a 2 = At, so that a = 2y/i, 
R-4 = —12/t = —6t and tdfo/dt = 1/2. Finally, the generating function for arbitrary 
Eulerian planar graphs, with a weight t per edge, and g2% per 2i-valent vertex and with a 
marked oriented edge (with a net effect of 2td/dt) reads 

E(t) = 2t% - 2t^ = -( - i?_ 4 - a 2 - 2t) 
w dt dt 2t v ; 

1 h M ( A - 12 ) 
{-j ~ J2 ^h-ri - a 2 - 2t) 



2t v t 

i=i 
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Let us finally introduce the function U(t) — a 2 /A. With this definition, eqn. ( |A.9| ) reduces 
to 

T-E^rV" 1 (A ' 13) 

while ( |A.12|) gives 

i=2 v 7 



i+1 



It is also interesting to compute the following derivative of E: 



|(tBW) = £ - 1 (A.15) 



(A. 14) 



where we have used ( |A.13| ) to eliminate dU/dt. This relation, due to Bender and Canfield 
33(1 , was given a combinatorial interpretation in Ref. 0, by establishing a bijection be- 



tween suitably marked Eulerian graphs and pairs of trees generated by the function U(t). 
Note that in this approach the Gaussian free energy /o is not subtracted from the free 
energy /, resulting in a shift E — > E + 1. 
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